The oxygen isotope effect in the ab-plane reflectance of underdoped YBa 2 Cu 3 07_j. 
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We have measured the effect of oxygen isotope substitution on the ab-plane reflectance of under- 
doped YBa2Cu307_s. The frequency shift of the transverse optic phonons due to the substitution 
of 16 by 18 O yields an isotope effect of the expected magnitude for copper-oxygen stretching modes 
with a = 0.5 ± 0.1. The reflectance shoulder at 400 - 500 cm -1 shows a much smaller exponent of 
a — 0.1 ± 0.1 in the normal state and a = 0.23 ± 0.1 in the superconducting state. These observa- 
tions suggest that the shoulder is of electronic origin and not due to a phonon mode as has been 
suggested recently. 

PACS numbers: 74.25.Kc, 74.25.Gz, 74.72.-h 



The absence of a substantial isotope effect on the tran- 
sition temperature of high temperature superconductors 
has been used to rule out phonons as the bosons respon- 
sible for coupling between the charge carriers Q. But 
this lack of T c shift with isotope substitution has failed 
to convince everyone. Many mechanisms in addition to 
the frequencies of the underlying phonon can affect T c . 
For example, among BCS superconductors, there are ex- 
amples where the isotope effect of T c is small and even 
positive due to the Coulomb interaction parameter fi* . 
whereas we know that the underlying boson frequencies, 
phonons in this case, would show a square root depen- 
dence on isotopic mass. Thus the isotope effect of T c is of 
marginal value in settling the question of phonon origin of 
superconductivity. A direct measurement of the isotope 
effect of the boson frequency is a much more convincing 
test. 

Recently, attention has been directed to the kink in 
the dispersion of the angle-resolved photoemission spec- 
tra (ARPES)@, |||4[ §. The simplest picture of the kink 
is that it is the result of the interaction of the charge car- 
riers with a bosonic excitation. This opens a new scatter- 
ing channel at the energy required to generate the excita- 
tion, and according to Kramers-Kronig relations applied 
to the complex self energy, a kink is formed in the energy 
dispersion In infrared spectroscopy the kink can be 
seen as an onset of absorption or a "knee" in low tem- 
perature reflectance spectra around 500 cm -1 , a result 
of enhanced scattering above this energy || 0]. The en- 
ergy of the knee does not vary among the various high 
temperature superconductors with the notable exception 
of the three-layer mercury compound where it is at a 
substantially higher energy, at 750 cm _1 [|). 

The knee has been associated with the superconduct- 



ing transition (in optimally doped YBa2Cu307_5 it can 
only be seen below T c ), although in most underdoped 
materials it is present in the normal state and has there- 
fore been related to the pseudogapQ. This has lead 
to the suggestion that the bosonic excitation responsi- 
ble for the knee is a possible candidate for the pairing 
excitation of high temperature superconductivity^, fL0| , 
similar to the role played by phonons in conventional 
superconductivity. Recently Lanzara et al. have pro- 
posed that certain zone boundary oxygen vibrations are 
responsible for the kink in ARPES dispersion curves^). 
A more conventional view is that the excitation is of mag- 
netic origin |l0j [H]. Suggestions that it is related to the 
sharp 41 meV neutron resonance have been challenged 
recently Q. 

An obvious experiment that will distinguish between 
these two models of the kink excitation is the oxygen iso- 
tope effect of the frequency of the kink. If the kink is due 
to phonons, its frequency should vary as the -0.5 power- 
law of the reduced mass of the mode. On the other hand, 
the frequency of a purely magnetic mode should be in- 
dependent of isotopic mass. Any such experiment will 
be difficult because the mode giving rise to the kink is 
broad. This is illustrated in Fig. la where we use stan- 
dard transport theory |b| to calculate the conductivity 
and the reflectance in a Fermi liquid system coupled to a 
bosonic mode. We model the kink excitation as a copper- 
oxygen stretching mode phonon of 12 cm" 1 width cen- 
tered at 516 cm -1 , interacting with the charge carriers 
with a plasma frequency tu p = 1.9 eV and a coupling 
constant A = 2.2. Fig. la shows that within this model, 
the complete replacement of 16 by 18 in the copper 
oxygen planes, results in the displacement of the knee in 
reflectance to lower frequency curve by 24 cm -1 , or al- 
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FIG. 1: Top curves (a), model calculation of the reflectance 
change on oxygen isotope substitution at 25 and 80 K for an 
copper-oxygen stretching mode phonon at a frequency of 516 
cm -1 coupled to the charge carriers with A = 2.2 and lo p = 1.9 
eV. Bottom curves (b) straight lines fit to reflectance experi- 
ments on isotopically separated YBa2Cu307_«. For each tem- 
perature, two curves are shown, the top curve (dashed) is for 
16 0, the bottom one (solid) for ls O. 



ternately, a reduction of reflectance in the spectral region 
above the knee by 0.7 %. We express a partial isotope ef- 
fect of the frequency of the mode v in terms of the isotope 
coefficient a defined in dvjv — —adm/m where m is the 
reduced mass of the copper-oxygen stretching vibration. 

To our knowledge, there are no published experiments 
on the isotope dependence of the frequency scale associ- 
ated with the pseudogap state or the kink in dispersion 
of the charge carriers. There are however several reports, 
based on NMR and NQR measurements, on the isotope 
dependence of the temperature T* associated with the 
onset of the pseudogap in YBa2Cu40s, a material that 
is naturally underdoped with a T c — 80 K. A zero iso- 
tope coefficient was reported by Williams et al. [ [Til , 
whereas Raffa et al^H^ found a value of a = 0.061 for T*. 
The same value, within experimental error, of a — 0.056 
was found for the superconducting transition tempera- 
ture T c , defined as a in dT c /T c — —adm/m where m 
is the mass of the oxygen isotope. Giant isotope effects 
have been reported on the relaxation rate of crystal field 
excitations JTg[ . 

For the experiment we used a large single crystal of 
YBa2Cu30*7_,5 which was broken in two pieces which 
were parallel annealed, one in 16 O and the other in 18 O 
to yield two underdoped samples. The weight changes 
of the crystals were consistent with a near complete iso- 
tope exchange. Neutron diffraction experiments showed 
that the occupancy of the chain sites was identical for the 
two samples. More importantly, the measurements also 



FIG. 2: The reflectance of YBa2Cu 3 O7-S, top curve, and 
isotopically substitute YBa2Cu3 8 07_s, lower curve. The fre- 
quency shift of two transverse optic phonons, Pi and P2 can 
be seen whereas two electronic excitations, Bi and B2 do not 
shift. For clarity the ls O curve has been shifted down by 0.02. 

show that the z-parameter of apical oxygen only shifts by 
-0.00003(4) c-axis units which is the opposite sign from 
what one would expect if the shift in T c was solely due 
to a difference in oxygen occupation (+0.00006) |l7). The 
uncertainty in the apical oxygen position gives an uncer- 
tainty in T c of 0.45 K. The change in shoulder frequency 
with Tc is approximately 2.5 cm /K giving an uncer- 
tainty in the shoulder frequency due to occupancy of the 
chain site of only 1.1 cm -1 . This is too small to compen- 
sate for a possible isotope shift expected due to phonons 
(see Fig. 1) of 24 cm -1 . 

The superconducting transition temperatures were T c 
= 67.6 K for 16 and T c = 66.7 K for 18 which 
yields a (partial) isotope coefficient a = 0.1. For a 
pure electron-phonon interaction in a monatomic system, 
a = 0.5. A small value of a is typical of high temper- 
ature superconductors jjj . The advantage of an under- 
doped sample is that it allows us to measure the kink 
frequency in both the normal and the superconducting 
state. The samples were twinned. 

The reflectance measurements were performed on 
freshly cleaved faces normal to the c-axis. To correct 
for irregularities of the surface, resulting from the cleav- 
age process, the sample was coated with a gold layer, 
evaporated in situfH^. 

All the measurements were carried out in a cold-finger 
flow cryostat with a minimum sample temperature of 25 
K. A Bruker IFS 66 v/S spectrometer was used in the far- 
infrared with He-cooled bolometer detectors below 700 
cm" 1 and an MCT detector up to 8000 cm -1 . The ab- 
solute reflectance measurements have an error of ± 0.5 
%. 



3 



Figure 2 shows the reflectance of the two crystals at 25 
K. Since the spectra are nearly identical on the scale of 
the figure, the 18 curve has been displaced downward by 
0.02 for clarity. We have normalized the two spectra in 
the 250 to 420 cm^ 1 region, below the reflectance shoul- 
der frequency, by subtracting 0.3 % from the reflectance 
of the 16 curve at 25 K and 0.75 % at 80 K. 

In Fig. 2 the reflectance minima marked Pi (bond 
bending) and P2 (bond stretching) are due to transverse 
optic phonons ||li|. The effect of the isotopic substitu- 
tion on their resonant frequencies can be seen clearly. 
From the isotopic masses one expects a frequency ratio 
(FR) Vis/^ie — 0.9428 for pure copper-oxygen stretch- 
ing vibrations. The observed values of FR at 25 K are 
0.956 and 0.953 for Pi and P2 respectively. The ex- 
pected values from the shell model of Chaplot et a/.]20|] 
are 0.954 and 0.9456 respectively. The zone boundary 
half breathing mode is predicted to have an isotope ra- 
tio of 0.9437, close to the pure oxygen vibration. Note, 
however, that the anomalously low frequency of the half- 
breathing mode is not correctly described by the model. 
The good agreement of FR with the theoretical values 
is consistent with complete in-plane isotopic substitution 
since the isotope effect on phonon frequencies is reduced 
for isotopic mixtures in a linear fashion plj]. 

The modes marked Bi and B2 do not shift significantly 
with isotopic substitution. Their frequencies and spectral 
weights suggest that they may be related to the modes 
discussed recently by Bernhard et al. [^2| . These authors 
found two modes in this spectral region with spectral 
weights over an order of magnitude larger than what is 
expected for phonons, and suggested that they may be 
of electronic origin. The absence of an oxygen isotope 
shift supports this picture although, it should be pointed 
out, low frequency modes generally involve heavier atoms 
in the unit cell and would not be expected to exhibit a 
complete oxygen isotope shift. While the modes Bi and 
B2 have been seen by previous workers [p2|, they were not 
seen in ultra-pure crystals)^ and arc likely activated by 
impurities. 

We obtained the frequency shift of the shoulder 
with isotope substitution in several ways. The most 
straightforward method would be to evaluate the second 
derivative of the reflectance which is closely related to 

a 2 F(f2)(H], |§ H S which has a P eak at the fre- 
quency of the mode responsible for the shoulder in the 
conventional Fermi liquid picture. We find this peak at 
520 cm -1 in the superconducting state and at 420 cm -1 
at 80 K, in the normal state. However, as a recent study 
shows, even with large crystals of Bi 2 Sr 2 CaCu208 which 
have very flat surfaces, the error of this method is too 
large to resolve the isotope effect (2^]. The source of this 
error is the well known growth of error with successive 
differentiation. To find a peak position in the a 2 F(il) 
spectrum, we must in effect evaluate the zero-crossing of 
the third derivative of the reflectance. 



To achieve the required resolution, we adopted the sim- 
pler method of evaluating the shift along the frequency 
axis of the slope portion of the reflectance in the 480 to 
620 cm -1 region. As Fig. 1 shows, within the Fermi liq- 
uid model, the isotope effect gives rise to a uniform shift 
of the curves along the frequency axis. We evaluated this 
shift in three different ways. The simplest was to fit a 
straight line to the reflectance in the slope region and 
to evaluate the frequency shift of this line for the two 
isotopes. This method gave a = 0.23 ±0.1 in the super- 
conducting state at 25 K and a = 0.1 ±0.1 in the normal, 
pseudogap state, at 80 K. To minimize uncertainties in 
the values of the absolute reflectance, all the curves were 
normalized in the 250 to 420 cm -1 region. In the second 
method, we again normalized the curves in the flat region 
and and evaluated the average shift in reflectance in the 
slope region and divided this by the slope of reflectance, 
dR/dcu to get the frequency shift. Within an experimen- 
tal error of a of 0.1, the results were the same as those 
obtained by the first method. 

We were concerned with the influence of the optic 
phonons Pi and P2 to our analysis. The third, more 
elaborate method, was designed to subtract the optic 
phonons from the spectra. We started by fitting the nor- 
malized reflectance to a series of six Lorentzian oscilla- 
tors that described the optical conductivity in the 250 to 
1000 cm -1 region. These included a Drude oscillator and 
the two phonons. We then recalculated the reflectance, 
without the phonons, and fitted a straight line in the 
slope region to the calculated model reflectance. The re- 
sulting isotope shift was identical to the one found by a 
direct fit to the experimental reflectance that included 
the phonons. 

The increased isotope dependence of the shoulder in 
the superconducting state may be related to the isotope 
effect of the superconducting transition temperature [Q. 
In our samples a(T c ) — 0.1±0.1. In several models of the 
shoulder at 520 cm -1 , its frequency in the superconduct- 
ing state is a sum of the normal state mode frequency 
and the superconducting gap frequency [|l0| If there 
is a proportionality between the superconducting transi- 
tion temperature and the gap frequency, we would expect 
an additional isotope effect of 0.1 in the superconducting 
state within the scope of these models. This is roughly 
what we observed. 

Another experiment that is difficult to reconcile with 
the phonon model for the shoulder is the infrared re- 
flectance of the three-layer mercury-based superconduc- 
tor which has a T c of 130 K. In this material the shoul- 
der frequency is at 750 cm -1 , some 45 % higher than 
in materials with a maximum T c in the 94 K region^. 
There are no direct measurements of the oxygen stretch 
frequency in this material. However, one can estimates 
this frequency from the general rule lu 2 oc a~ 3 where a 
is the lattice parameter p9j |30| |. The lattice parameter 
of the three-layer mercury compound is 3.854 at optimal 
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doping, as compared to 3.827 in YBa 2 Cu307_5. The 
calculated phonon frequency shift, based on the lattice 
parameter change, is -0.5 %. This is to be compared to 
the observed shift of the shoulder of +45 %. 

In summary, we have measured the effect of substi- 
tuting the oxygen isotope 18 O for 16 on the infrared 
reflectance spectrum of underdoped YE^C^C^-a. We 
find that two peaks, previously assigned to phonons, 
show the expected isotope shift with a coefficient close 
to the theoretical value of 0.5, whereas the reflectance 
sho ulder has a much smaller coefficient of 0.1 ±0.1 in the 
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nor mal state which increases to 0.23 ± . 1 in the super - 



conducting state. Thus, within experimental error, we 
can rule out any simple phonon model involving copper- 
oxygen stretching modes as an explanation for the re- 
flectance shoulder in the normal state. The small isotope 
shift of the shoulder suggests that it is largely of elec- 
tronic origin. 

The work at McMaster University and the University 
of Alberta was supported by the Canadian National Sci- 
ence and Engineering Research Council. We thank J. P. 
Carbotte, D.N. Basov, A. Lanzara. M.R. Norman and 
Z.-X. Shen for useful discussions. 



* Present address: Institute of Physics, Chinese Academy 
of Sciences, P.O.Box 2711, Beijing 100080, P. R. China 

t The Canadian Institute of Advance d Research; Elec- 
tronic address: timusk@mcmaster.ca 

1 Laboratoire Leon Brillouin, CE Saclay, F-91191 Gif-sur- 
Yvette, France; Present address: II. Phys. Inst., Univ. zu 
Koln, Zuelpicher Str. 77 D-50937 Koln, Germany 

§ Present address: Walther Meissner Institut, Bayerische 
Akademie der Wissenschaften, D-85748 Garching, Ger- 
many 

[1] For a review, see J. P. Franck in Physical Properties of 

High Temperature Superconductors III, D.M.Ginsberg, 

editor, (World Scientific, Singapore, 1994) pp. 189-293. 
[2] A. Kaminski, J. Mesot, H. Fretwell, J.C. Campuzano, 

M.R. Norman, M. Randeria, H. Ding, T. Sato, T. Taka- 

hashi, T. Mochiku, K. Kadowaki, and H. Hoechst, Phys. 

Rev. Lett. 84, 1788 (2000). 
[3] T. Valla, A.V. Fedorov, P.D. Johnson, Q. Li, G.D. Gu, 

and N. Koshizuka, Phys. Rev. Lett. 85, 828 (2000). 
[4] P.V. Bogdanov, A. Lanzara, S.A. Kellar, X.J. Zhou, E.D. 

Lu, W.J. Zheng, G. Gu, J. -I. Shimoyama, K. Kishio, H. 

Ikeda, R. Yoshizaki, Z. Hussain, and Z.X. Shen, Phys. 

Rev. Lett. 85, 2581 (2000). 
[5] P.D. Johnson, T. Valla, A.V. Fedorov, Z. Yusuof, 

B.O. Wells, Q. Li, AR. Moodenbaugh, G.D. Gu, N. 

Koshizuka, C. Kendziora, Sha Jian, and D.G. Hinks, 

Phys. Rev. Lett. 87, 177007 (2001). 
[6] D.N. Basov, R. Liang, B. Dabrowski, DA. Bonn, 

W.N. Hardy,and T. Timusk, Phys. Rev. Lett. 77, 4090 

(1996). 



[13; 
[14; 
[is; 

[ie; 

[17. 

[is; 

[19 

[20 
[21 

[22 

[23 

[24 
[25 

[26 

[27 

[28 

[29 

[30 



A.V. Puchkov, D.N. Basov, and T. Timusk, J. Physics, 

Condensed Matter, 8, 10049-84, (1996). 

J.J. McGuire, M. Windt, T. Startseva, T. Timusk, D. 

Colson, V. Viallet-Guillen, Phys. Rev. B 62, 8711 (2000). 

A. Lanzara et al. Nature 412, 510 (2001). 

J. P. Carbotte, E. Schachinger, and D.N. Basov, Nature 

401, 354 (1999). 

P. Bourges, in The gap Symmetry and Fluctuations 
in High Temperature Superconductors, ed. J. Bok, 
G. Deutscher, D. Parrina. and S. A. Wolf, Plenum Press 
1998, preprint |cond-mat/9901333; . 

Hae- Young K ee, Steven A. Kivelson, G. Aeppli, cond 



mat /0 110378 



S.V. Shulga, O.V. Dolgov, and E.G. Maksimov, Physica 
C 178, 266 (1991). 

G.V.M. Williams, J.L. Tallon, J.W. Quilty, H.J. Trodahl, 

and N.E. Flower, Phys. Rev. Lett. 80, 377 (1998). 

F. Raffa, T. Ohno, M. Mali, J. Roos, D. Brinkmann, 

K. Conder, and M. Eremin, Phys. Rev. Lett. 81, 5912 

(1998). 

D. Rubio Temprano, J. Mesot, S. Janssen, A. Furrer, K. 
Conder, and H. Mutka, Phys. Rev. Lett. 84, 1990 (2000). 
The structure of the two crystals was studied on the four 
circle diffractometer D9 at the Institut Laue Langevin 
(Grenoble), where data sets of 560 independent Bragg 
reflection intensities were collected. Details of this study 
will be given elsewhere. 

C.C. Homes, MA. Reedyk, D.A. Crandles, and T. 
Timusk, Applied Optics 32, 2976 (1993). 

C. C. Homes, A. W. McConnell, BP. Clayman, D.A. 
Bonn, Ruixing Liang, W.N. Hardy, M. Inoue, H. Negishi, 
P. Fournier, and R. L. Greene, Phys. Rev. Lett. 84, 5391 
(2000). 

S.L. Chaplot, W. Reichardt, L. Pintschovius, and 
N. Pyka, Phys. Rev. B 52, 7230 (1995). 

D. J. Montgomery and JR. Hardy, in Lattice Dynam- 
ics, Proc. of the Int. Conf. held in Copenhagen 1963, 
R.F. Wallis editor, Pergamon Press, Oxford (1965) p. 
491 

C. Bernhard, T. Holden, J. Humlfcek, D. Munzar, A. 
Golnik, M. Klaser, Th. Wolf. L. Carr. C. Homes. B. 



Keimer, and M. Cardona, sond-mat /010757 



C. C. Homes, D. A. Bonn, Ruixing Liang, W. N. Hardy, 

D. N. Basov, T. Timusk, and BP. Clayman, Phys. Rev. 
B 60, 9782 (1999). 

PB. Allen, Phys. Rev B 3, 305 (1971). 

B. Farnworth and T. Timusk, Phys. Rev B 10, 2799, 

(1974). 

F. Marsiglio, T. Startseva, J. P. Carbotte, Phys. Lett. A 
245, 354 (1998). 

Ar. Abanov, A.V. Cubukov, and J. Schmalian, Phys. Rev 
B 63, 180510 (2001). 

J.J. Tu, C.C. Homes, G.D. G u, D.N. Basov, S. M. 



Loureiro, R.J. Cava, M. Strongin, cond-mat/0104208 
S. Tajima, T. Ido, S. Ishibashi, T. Itoh, H.Eisaki, 
Y. Mizuo, T. Arima, H. Takagi, S. Uchida, Phys. Rev 
B 43, 10496 (1991). 

L. Pintschovius, J.M. Bassat, P. Odier, F. Gervais, 
G. Chevrier, W. Reichardt, and F. Gompf, Phys. Rev. 
B 40, 2229 (1989). 



